Quantitative analysis of single particle trajectories: mean maximal excursion method 
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An increasing number of experimental studies employ single particle tracking to probe the 
physical environment in complex systems. We here propose and discuss new methods to analyze 
the time series of the particle traces, in particular, for subdiffusion phenomena. We discuss the 
statistical properties of mean maximal excursions, i.e., the maximal distance covered by a test 
particle up to time t. Compared to traditional methods focusing on the mean squared displacement 
we show that the mean maximal excursion analysis performs better in the determination of the 
anomalous diffusion exponent. We also demonstrate that combination of regular moments with 
moments of the mean maximal excursion method provides additional criteria to determine the 
exact physical nature of the underlying stochastic subdiffusion processes. We put the methods to 
test using experimental data as well as simulated time series from different models for normal and 
anomalous dynamics, such as diffusion on fractals, continuous time random walks, and fractional 
Brownian motion. 
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I. INTRODUCTION 

The history of stochastic motion may be traced back 
to the writings of Titus Lucretius, describing the battling 
of dust particles in air Later, irregular motion of sin- 
gle coal dust particles was described by Jan Ingenhousz 
in 1785 Robert Brown in 1827 reported the jittery 
motion of small particles within the vacuoles of pollen 
grains Q. Possibly the first systematic recording of ac- 
tual trajectories was published by Jean Perrin, observing 
individual, small granules in uniform gamboge emulsions 
0. Yet apparently the first experimental study based 
on the time series analysis of single particle trajectories 
is due to Nordlund who tracked small mercury spheres 
in water Today single trajectory analysis is a com- 
mon method to probe the motion of particles, notably, 
in complex biological environments |6l-ll5| . 

Typically a diffusion process in d dimensions is charac- 
terized by the ensemble averaged mean squared displace- 
ment (MSD) 
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Here we assumed spherical symmetry and an isotropic 
environment, such that P{r,t) is the probability density 
to find the particle a (radial) distance r away from the 
origin at time t after release of the particle at r = at 
time t = 0. In equation ([1} we introduced the anomalous 
diffusion exponent a. In the limit a = 1 we encounter 
regular Brownian diffusion. For other values of a the 
associated diffusion is anomalous: the case < a < 1 
is called subdiffusion while for a > 1 the process is su- 
perdiffusive [l^. In this work we focus on subdiffusivc 



processes. In equation ([T]) the generalized diffusion co- 
efficient is of dimension [Ka\ = cm^/sec". Subdiffusion 
of the form ([l} is found in a variety of systems, such 
as amorphous semiconductors [l7j . tracer dispersion in 
subsurface acquifcrs [l^, or in turbulent systems [T9| . 

In fact, subdiffusion was found from observation of sin- 
gle trajectories in a number of biologically relevant sys- 
tems: For instance, it was shown that adeno-associated 
viruses of radius « 15 nm in a cell perform subdiffusion 
with a = 0.5 . . . 0.9 Q. Fluorescently labeled messenger 
RNA chains of 3000 bases length and effective diameter 
of some 50nm subdiffuse with a w 0.75 0. Lipid gran- 
ules of typical size of few hundred nm exhibit subdiffusion 
with a w 0.75 . . . 0.85 [i-[ll|; and the diffusion of telom- 
eres in the nucleus of mammalian cells shows a « 0.3 at 
shorter times, and a w 0.5 at intermediate times [12[. A 
study assuming normal diffusion for the analysis of track- 
ing data of single cell nuclear organelles shows extreme 
fluctuations of the diffusivity as function of time along 
individual trajectories, possibly pointing to subdiffusion 
effects [l3| . In vitro, subdiffusion was measured in pro- 
tein solutions [l^ and in reconstituted actin networks 
. Molecular crowding is often suspected as a cause of 
subdiffusion in living cells [lO, [U ■ 

Currently one of the important open questions is what 
physical mechanism causes the subdiffusion in biological 
systems. Single particle tracking is expected to provide 
essential clues to answer this question. Thus, recently 
a method has been suggested based on the statistics of 
first passage times, i.e., the distribution of times it takes 
a random walker to first reach a given distance from its 
starting point. This quantity has been shown to be a 
powerful tool to discriminate between CTRW and dif- 



2 



fusion on fractals [22|, |23|. However such an analysis 
requires a huge amount of data to be statistically rel- 
evant. Fluorescence correlation spectroscopy (FCS) has 
also been proposed to identify the physical mechanism of 
subdiffusion [2J]; but this approach is based on an indi- 
rect observable, the fluorescence correlator, which is not 
directly comparable with analytical results; moreover this 
method needs to fit three parameters to a single curve. 
We here present a new method, that is based on analyti- 
cal results. Our approach is demonstrated to enable one 
to extract more, and more accurate, information from a 
set of single particle trajectories. 

A typical single particle tracking experiment provides 
a time series r(t) of the particle position from which one 
may calculate the time averaged mean squared displace- 
ment 
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Here T denotes the overall measurement time, and A 
is a lag time defining a window swept over the time se- 
ries. For a Brownian random walk with typical width 
((5r^) of the step length and characteristic waiting time 
T between successive steps, we recover the time average 
(s^{A,T)^ = 2dKiA, where the diffusion constant be- 
comes Ki = ((5r^)/[2dT]. In this case the time average 
provides exactly the same information as the ensemble 
average. Note that this is not always the case when the 
dynamics is subdiffusive [25l - l27j . 

Using time averages to analyze the behavior of a single 
particle is an elegant method, in particular, to avoid er- 
rors from averages over particles with nonidentical phys- 
ical properties. However in many cases the actual trajec- 
tories arc too short to allow one to extract meaningful 
information from the time average. Moreover, in cases 
where the subdiffusion is governed by a CTRW with di- 
verging characteristic waiting time the values of the mo- 
ments. and therefore their ratios, become random quan- 
tities [2^ [2^ . Using the ensemble average prevents this 
problem. We therefore consider herein ensemble averages 
calculated directly from measured trajectories. In par- 
ticular we present an analysis based on a mean maximal 
excursion statistics. It will be shown that this method 
provides relevant information on the system, complemen- 
tary to results from analysis of regular moments. More- 
over we demonstrate that the mean maximal excursion 
method may obtain more accurate information about the 
dynamics than the typically measured mean squared dis- 
placement ([T|). 

In what follows we present the theoretical background 
of the mean maximal excursion analysis and discuss how 
different dynamic processes can be discriminated. We 
then discuss how to apply these methods in practice, in- 
cluding the analysis of some recent single particle track- 
ing data. 



II. MATERIALS AND METHODS 

As a benchmark for our quantitative analysis we here 
define the three most prominent approaches to subdiffu- 
sion. Physically these processes are fundamentally differ- 
ent, while they all share the form ([1]) of the mean squared 
displacement. In the supplementary material we provide 
details on how we simulate the time series based on the 
stochastic models. 

(i) Continuous Time Random Walk (CTRW). CTRW 
defines a random walk process during which the walker 
rests a random waiting time, drawn from a probability 
distribution, between successive steps [17| . If the density 



of waiting times is of the long tailed form 
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for < a < 1, the mean waiting time t%p{t)dt di- 
verges, and the resulting process becomes subdiffusive 
with mean squared displacement ([T]). The exponent a 
from the waiting time density ([3]) is then the same as 
in equation ([!]). If the variance of the associated jump 
lengths is again ((5r^), the generalized diffusion coeffi- 
cient becomes = (5r^)/(2(ir). Waiting times with 
such power-law distribution were, for instance, observed 
for the motion of probes in a reconstituted actin network 
[l5| . CTRW is used in a wide variety of fields, ranging 
from charge carrier motion in amorphous semiconductors 
(TtI , over tracer diffusion in underground aquifers [l^ , up 
to weakly chaotic systems [l^. 

(ii) Diffusion on fractals. A random walker mov- 
ing on a geometric fractal, for instance, a percolation 
cluster near the percolation threshold, meets bottlenecks 
and dead ends on all scales, similar to the motion in a 
labyrinth. This results in an effective subdiffusion in the 
embedding space. While the fractal dimension df char- 
acterizes the geometry of the fractal, the diffusive dy- 
namics involves an additional critical exponent, the ran- 
dom walk exponent d^ {d^ > 2). The latter is related 
to the anomalous diffusion exponent through a = 2/d^ 
[28j . Fractals can be used to model complex networks, 
and have recently been suggested to mimic certain fea- 
tures of diffusion under conditions of molecular crowding 
j29l [30| . We will use for the theoretical descriptions the 
dynamical scheme of reference [sij . 

(Hi) Fractional Brownian Motion (FBM). FBM was 
introduced to take into account correlations in a random 
walk: the state of the system at time t is influenced by the 
state at time t' < t. In the FBM model this is achieved by 
passing from a Gaussian white noise dB{t) to fractional 
Gaussian noise 
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dB{T)], (4) 



where the Hurst exponent < i7 < 1 is connected to 
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the anomalous diffusion exponent by a = 2H . FBM 
therefore describes both subdifFusion and superdiffusion 
up to the baUistic hmit a = 2. FBM is used to describe 



the motion of a monomer in a polymer chain [32[ or single 
file diffusion [s^. FBM has recently been proposed to 
underlie the diffusion in a crowded environment |24l| . The 
autocorrelation function of FBM in ID reads 1341 
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and for ti = ^2 we recover the mean squared displacement 
([1]). Following reference [13, we extend FBM to several 
dimensions such that a d-dimensional FBM of exponent 
H \s a. process in which each of the coordinates follows 
a one-dimensional FBM of exponent H . The resulting 
d-dimensional FBM still satisfies ([T]), with a = 2H . 



III. RESULTS 

The parameters in the three simulation models are cho- 
sen to produce the same anomalous diffusion exponent 
a = 0.70. Using only the classical analysis based on the 
MSD ll]), one could not tell which model was used to 
create the data. We discuss here how additional observ- 
ables allow one to extract a more accurate value of this a 
exponent, and how they may be used to distinguish the 
microscopic stochastic mechanisms. 



A. Mean maximal excursion (MME) approach 

A power law fit to the classical MSD ([T]) provides the 
magnitude of the anomalous diffusion exponent a. We 
here show that the MME method is a better observable 
to determine a. The maximal excursion is the greatest 
distance r, that the random walker reaches until time t. 
This quantity is averaged over all trajectories, to obtain 
the MME second moment 
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where Pr (rmax — fo,t) is the probability that the maxi- 
mal distance from the origin that is reached up to time 
t, is equal to tq. The MME second moment ^ scales 
like t", as shown in reference [l^ for fractal media, and 
derived in the supplementary material for a CTRW pro- 
cess. 

For FBM this quantity is not known, similar to the 
first passage in other than a semi-infinite domain in ID. 
However, one can still use the MME method to numeri- 
cally analyze data created by an FBM process, as shown 
below. 

Why is the MME second moment better than the more 
standard MSD? The ratio 7 = ax{t)/{X{t)) of the stan- 
dard deviation axit) = y/{{X{t) - {X{t))y) versus the 
mean is a measure for the dispersion around the center 



of the distribution (first moment). A lower ratio means 
that the random variable has a smaller spread around 
its mean. This will produce a smoother average and 
thus a more accurate fit as the larger number of data 
points closer to the average value receive a higher rel- 
ative weight. Indeed, for regular Brownian motion the 
ratio is smaller for the MME second moment than for 
the regular second moment, the time independent val- 
ues being 7(MSD)/7(MME) = 1.61, 1.44, and 1.34 for 
one, two, and three dimensions. The MME method is 
therefore expected to non-negligibly outperform the MSD 
method. Details of this calculation are presented in the 
supplementary material. For diffusion on a fractal, the 
ratio 7(MSD)/7(MME) also grows with decreasing frac- 
tal dimension, being always greater than 1. For a CTRW 
the ratio 7(MSD)/7(MME) diminishes as well with de- 
creasing a, reaching its lowest value at a = 0. But it is 
always larger than 1 in dimensions d = 1, 2, 3. 

Another way to characterize the dispersion of the 
MME method versus regular moments is the ratio of the 
fourth moment versus the second moment of the respec- 
tive distribution: (i) For a random walk on a fractal. 



approximated by the dynamical scheme of reference [31 
the MME moments become 13(1 
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where the prefactor is given through 
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Here /„ is the modified Bessel function of the first kind. 
The regular moments satisfy an analogous relation [sij , 



V{a[k + df]l2) 
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The ratios (''max)/ (''max) ^ {r^)/{r'^y are therefore 
time independent numerical constants. Note that above 
expressions also contain the limiting case of Brownian 
motion (integer dimension, and a = 1). In the latter case 
the associated values are listed in table HI demonstrating 
again that the MME distribution is more concentrated 
and therefore more amenable to parameter extraction by 
fitting, see also the discussion below. 

(ii) For FBM, the regular moments are obtained from 
the Brownian ones by simple replacement of time thyt°'. 
Since the regular moment ratios arc time independent 
we find exactly the same values as in the Brownian case. 
The MME moments arc not known analytically, so we 
performed numerical simulations to get an estimate of 
these quantities. A surprising result is that the MME 
moments (j'max) are proportional to t''" but with a 
new exponent a' > a. 

We discuss these results in detail in the supplementary 
material, finding a linear correlation (R^ > 0.999 for 10 
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1.36 
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2.78 
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TABLE I: Ratios of fourth moment versus the square of 
the second moment for normal moment statistics and MME 
statistics. We hst normal Brownian motion {a — 1) and 
CTRW subdiffusion with a = 1/2. The MME distribution 
is narrower and therefore more amenable for data fitting in 
all cases. 

points) between the two exponents: 

a' « 0.156 ± 0.005 + (0.849 ± 0.008)a (10) 

We note that for Brownian motion (a = 1), we retrieve 
the classical result a' = a. We also obtained an expres- 
sion for the MME moment ratio, (?'max)/(''max)^j 2D 
(i?2 > 0.99 for 10 points): 

/r-4 \ /n-x i-42±0.0l 

ff^ ^ (1.05 ±0M)(^) +(1.10 ±0.01). 

(11) 

We note that solely focusing on the determination of a' 
from the second MME moment may lead to an overes- 
timation of the anomalous diffusion exponent if the mo- 
tion is governed by FBM and a' is not converted to a 
via relation ([TT]) . It is therefore important to also evalu- 
ate the complementary criteria such as the mean squared 
displacement and the moment ratios. 

(iii) In the case of CTRW subdiffusion we profit from 
the fact that in Laplace space we can transform the prob- 
ability density and the moments of normal Brownian 
motion into the corresponding CTRW subdiffusion so- 
lution by so-called subordination [l^, [s^]. In practice 
this means that we can replace s by Kis" / Ka where s is 
the Laplace variable conjugated to time t. We obtain the 
ratio for both regular moments and MME statistics from 
the Brownian result, however, with different pre-factors 

k r{k/2 + 1) , k \ 

(^max)cTRW = rXofcT^TT) 

Table U shows the results for a — 1/2. 

The moment ratios (''max)/('"max)^ ^iid {r^) / {r"^)"^ are 
useful observables. Once we determine the anomalous 
exponent a from fit to the MSD or the second MME mo- 
ment we can use the moment ratios to identify the pro- 
cess. If the moment ratio for a subdiffusion process with 
< a < 1 is the same as for Brownian motion we are 
dealing with an FBM process. If the value matches the 
one for CTRW subdiffusion for the given a we verify the 
CTRW mechanism. Finally, we can identify the remain- 
ing possibility, i.e., diffusion on a fractal: The obtained 



numerical value for the ratio allows us, in principle, to 
deduce the underlying fractal dimension d/, using the 
predicted values of equation ^ and dH). We will discuss 
below how reliable such classifications are. 



B. Determination of the fractal dimension df 

Finally we establish a criterion to distinguish diffusion 
on a fractal from CTRW and FBM subdiffusion. We 
know that the probability density for a diffu sing particle 
on a fractal satisfies the scaling relation [H, l39| 

P(r,i)=t-"'^//2p(^,l). (14) 

The same relation holds for a CTRW or a FBM if we 
replace df hy the Euclidian dimension. Let us focus on 
the probability to be in a growing sphere of radius rot"/^. 
Then 

Pr^r<^o^"/^^j = / r'^-'^ P{r,t)dr 

= A{ro)t"^'^-'^f^/'^. (15) 

Since the exponent a is known from the second MME 
moment fit we can extract df from above relation. 

C. Summary 

Collecting the results from this section we come up 
with the following recipe to analyze diffusion data ob- 
tained from experiment or simulation, compare also the 
results summarized in table |TT1 

(1) Obtain the anomalous diffusion exponent a from 
power law fit to MSD and second MME moment. Dif- 
ferent subdiffusion mechanisms can the be determined as 
follows: (2) Diffusion on a fractal has regular and MME 
moment ratios, that depend on both a and the fractal 
dimension df. The fractal dimension is smaller than the 
embedding Euclidean dimension. (3) CTRW subdiffu- 
sion has regular and MME moments that depend on the 
anomalous diffusion exponent a. The ratios are larger 
than the corresponding Brownian quantities. The prob- 
ability to be in a sphere growing like t"/^ is constant. 
(4) FBM has the same ratios for regular moments as 
Brownian motion. The MME second moment exponent 
is greater than a, and the MME ratio is smaller than the 
Brownian one. The probability to be in a sphere growing 
like t"/^ is constant. 



IV. DISCUSSION 

We now turn to the question how experimental data 
can be analyzed by help of the tools established above. 
In a typical experiment a small particle is tracked by a 
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Second moment (regular, MME) 


Ratio (regular, MME) 


Growing spheres 


BM 


(oc t,oc t) 


(2,1.49), eq. © and ® 


Pr(r <ro^"/^^) = 


Fractals 


(oc t",oct") 


(< 2, < 1.49), eq. (0 and Q 


Pr(r <ror/^^) oct'^(2-''/'/2 


CTRW 


(oc r,ocr) 


(> 2, > 1.49), eq. ((T2| and dTSj 


Pr (r < rot"/2,t) = Tlo 


FBM 


(oc t",oc t"'), eq. ^ 


(2, < 1.49), eq. (0 


Pr (r <ror/2,t) = Ao 



TABLE II; Test for 2D trajectories in a free environment, and equation references for other dimensions. For each model, 
Brownian motion (BM), diffusion on fractal, continuous time random walk (CTRW), and fractional Brownian motion (FBM), 
the second column lists the scaling behavior of the second regular and MME moments ((r^) and {r'^-^^^)); the third column 
shows the relative values of the regular and MME ratio ({X'^) / (X^)'^); and the fourth column contains the scaling laws of the 
probability, at time t, to be in a sphere growing like 1°''''^. 



microscope, the motion being projected onto the focal 
plane (2D), to produce a time series r(t) = {x{t),y{t)) 
of the particle positions. Given a set of N trajectories 
ri(t), with Hi steps in trajectory i, we first calculate the 
distances to the starting point, 

n{t) = ^ [x,{t) - x,iO)f + [y.,{t)-y,m\ (16) 

in the 2D projection of the motion monitored in the ex- 
periment. The propagator is not directly accessible in an 
experiment. However, division of the number of trajecto- 
ries being at r for a given time t in the 2D projection, by 
the total number of trajectories of length m > t, leads to 
a good estimate of P(r, t) . We can therefore transform all 
the previous integrals defining the moments into discrete 
sums, and apply above methods. 
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A. Regular and MME moments 

In discrete form the fcth order moments become 

AA(t) 

<^'W)-AmyE-'W (17) 

and 

1 r 

('^LxW)-^E(o^^,^^,{'^^(0}j , (18) 

for regular and MME statistics, respectively. Here J\f{t) 
is the number of trajectories that are at least t steps long. 

Note that the discrete MME moments defined here do 
not correspond exactly to the theoretical definition pro- 
vided before. In fact, we do not have access to the whole 
trajectory, but only some sample points of it, with a given 
time step between two consecutive frames. The real Tmax 
may be reached in between two frames, and therefore 
would not be observed. However after sufficiently long 
time the difference between the discrete estimate calcu- 
lated here and the real value from the continuous trajec- 
tory becomes sufficiently small. 

Figure [T] shows the result of fits of the MSD and the 
second MME moment to simulated data according to the 



FIG. 1: MSD {r^(t)) and second MME moment (r,^,ax) as 
function of time t (arbitrary units) for the three simulated 
time series (1,000 trajectories of 100 steps each), each with 
anomalous diffusion exponent a = 0.7. The power law fits 
produce: (i) for 2D percolation data a — 0.64 (MSD, depicted 
by black x) and a = 0.73 (MME, black A); (ii) for CTRW 
data a = 0.67 (MSD, red x) and a = 0.71 (MME, red A), 
(iii) for FBM data a = 0.72 (MSD, green x) and a' = 0.79 
(MME, green A, expected value a' f» 0.74). 



three subdiffusion models, all with anomalous diffusion 
exponent a = 0.7. Indeed the MME method performs 
somewhat better. We should note that these simulation 
results are fairly smooth, and therefore we would not ex- 
pect a significant difference between the two methods, 
in contrast to the results on the experimental data be- 
low. Also note that we chose different anomalous diffu- 
sion constants Ka to be able to distinguish the different 
curves in figure [T] Of course, this does not influence the 
quality of the fit of the anomalous diffusion exponent a. 

Let us now turn to the moment ratios {t'^) / {r'^)'^ and 
(^max)/(''max)^- mentioned above some care has to be 
taken with the latter: only the long time values have a 
physical meaning. In fact, for the first frame, the mo- 
ment estimate ('r^ax) is exactly (r^), because of the dis- 
crete time step. After few dozens of frames, the estimate 
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(''max) converges toward its correct value, and the ratios 
become meaningful. 

In figure [2] we show a plot of the moment ratios. The 
convergence to a constant value attained at sufficiently 
long times is distinct. The ratios are those predicted 
for both CTRW and FBM, where the simulation is per- 
formed in a free environment. For diffusion on a perco- 
lation cluster, we observe a deviation from the predic- 
tion, due to the confinement of the diffusion for this set: 
the propagator does not converge toward the free space 
propagator, but toward the stationary distribution. We 
note that these ratios are clearly distinguishable between 
regular and MME moments, but also between the three 
simulations sets. Knowing the a value from the previous 
power law fit of MSD or second MME moment, those 
ratios are already a good indication of the underlying 
stochastic process. Since the difference between CTRW 
and diffusion on a fractal is not too large, we use the 
method of a growing sphere to see whether we can dis- 
criminate more clearly between those two mechanisms. 

B. Growing sphere analysis 

Let us turn to the probability to find the particle at 
time t in a (growing) sphere of radius rot"^^. Here ro is 
a free parameter. It should be chosen sufficiently large, 
such that for a given trajectory the probability to be 
within the sphere is appreciably large. At the same time 
it should not be too large, otherwise the probability to 
be within the sphere is almost one. Choosing a small 
multiple of (r(t = 1)) appears to be a good compromise. 
The probability to be inside the sphere then becomes 

Pr (r < ro^/^) ^ {r.,{t) - r.^'^) . (19) 

^ ' i—l 

Here 9(r) is the Heaviside function, that equals 1 if r > 0, 
and if r < 0. We expect the scaling oc To fit 

the fractal dimension df we need the anomalous diffusion 
exponent a as input. We used the value extracted from 
the second MME moment fits. The direct plot of the 
probability is quite easy to interpret: if the probability is 
constant, then d = df, if it grows slowly, then d> df, and 
the support is fractal {df ^ d). The dimension d here is 
the dimension of the trajectories (c? = 2 in our examples 
due to the projection onto the focal plain). In figure [31 
we see clearly that for CTRW and FBM the probability 
is approximately constant, and that for the diffusion on 
a percolation cluster, it grows with time, indicating that 
df < d, as \i should be. 

C. Experimental data 

We analyse experimental single particle tracking data 
showing that such time series are sufficiently large to ap- 
ply the analysis tools developed herein. 
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FIG. 2: Regular and MME moment ratios (r'')/(r^)^ and 
{''max)/(''max)^ as functlon of time (a.u.) for the three sim- 
ulated sets (diffusion on a fractal, FBM, and CTRW). Each 
set consists of 1,000 trajectories with 100 steps each. 
Black A: MME ratio for the diffusion on a 2D percolation 
cluster; the data do not converge to the expected value 1.29 
(black horizontal line). The same behavior is observed for the 
regular moment ratio (black -I-), for which the expected value 
is 1.77 (short black line). This discrepancy is likely due to 
the confinement of the percolation cluster on a 250 x 250 net- 
work: the random walker quickly reaches the boundaries, and 
the convergence occurs toward the equilibrium distribution, 
not toward the free space propagator. 

Red A: MME ratio for the CTRW process, converging to 1.97 
(red horizontal line). We also plot the regular moment ratio 
(red +); these are more irregular and converge to 2.66 (short 
red line). 

For FBM, the MME ratio (green A) converges to the esti- 
mated value of equation Hll|) . 1.33 (green horizontal line), 
and the regular ratio (green +) oscillates around the Brown- 
ian value 2 (short green line). 



The first data set (see supplementary material) con- 
tains 67 trajectories with up to 210 steps length of quan- 
tum dots diffusing freely in a solvent. The expected be- 
havior is regular Brownian motion. The data set is quite 
small and we show that MME moments are better ob- 
servables than regular moments. We plot the MSD as a 
function of time in figure 21 and fit the data by a power- 
law cx i". This fit provides an anomalous diffusion co- 
efficient a = 0.81. The fit based on the second MME 
moment returns the value a = 1.02, an almost perfect 
reproduction of the expected value a — \. The much 
better result of the MME method is due to the lower 
dispersion around the mean of the MME statistics, as 
discussed in the supplementary material. In figure HI it 
can be appreciated that the large outlier in the MSD 
statistics at around t = 0.7 sec is responsible for the low 
a value. At longer times also the MSD follows normal 
diffusion. This analysis demonstrates that the MSD in 
this case would lead to a large deviation from the cx- 
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FIG. 3: Probability to be in a growing sphere of radius rat"^'^ 
as function of t"''^ for the three simulated sets (a.u.). This 
analysis is based on the previously fitted values of a. Results: 
(i) 2D critical percolation (black x) produces d — d/~0.11, 
i.e., df ^ 1.89 (exact value 91/48 f» 1.896). (ii) The CTRW 
set (red x) gives d — df ~ 0.01 instead of 0, and the FBM set 
(green x) leads to d — d/ ~ —0.004 instead of 0. 
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FIG. 4: Analysis of an experimental set of 67 trajectories, 
the longest consisting of 210 points, for quantum dots freely 
diffusing in a solvent. MSD (black x), fitted by a power 
law with exponents a — 0.81 (red line). We also show a fit 
with fixed exponent a = 1 (green line, expected behavior for 
Brownian motion). MME (blue x), fitted by a power law 
(red line, a = 1.02). Time is in seconds, distances are in 
Hm^. Inset: double-logarithmic plot of the same data. 



pected value, and thus to the erroneous eonclusion that 
the observed motion were subdifFusive, while the MME 
analysis perforins much more reliably. 

The second set of data was obtained from video track- 
ing of 8 different lipid granules moving in yeast cells. 
Since we had few long trajectories, before an ensemble 



FIG. 5: Lipid granules diffusing in a yeast cell: 8 trajecto- 
ries, between 5,515 and 19,393 frames long. Log-log plot of 
the time-averaged MSD as a function of lag time (continuous 
lines), and Aot"'* (dotted lines). Time is scaled in seconds, 
and the time averaged MSD in fj,m^ . 



average, we first directly analyzed the 8 trajectories us- 
ing the time-averaged MSD We obtain a distinct 
subdiffusive behavior with an exponent close to 0.4, as 
demonstrated in figure [S] Each trajectory corresponds to 
a given granule. It is interesting to see that the data ex- 
hibit a scatter in amplitude and considerable local varia- 
tion of slope. Such features were also observed previously, 
see, for instance, references 0, [l3|- They may possibly 
be related to ageing effects Ho]. We also note that one 
of the curves shows a much steeper slope than the oth- 
ers. We extended the time-average analysis to the second 
MME moment 
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and again obtained a clear subdifFusive behavior, but 
with an exponent close to 0.5, as demonstrated in fig- 
ure ini Once again, we have a scatter in amplitude. The 
initial slope variation {0 < t < 10) is due to the inac- 
curacy in the MME estimation when there are only few 
frames to average. A greater exponent for MME than 
for regular moment could be due to an inaccuracy in the 
fit. However, it may indeed point toward an underlying 
FBM process. 

In order to gain more insight into the diffusion mecha- 
nism producing this subdiffusion behavior, we applied the 
methodology detailed above. Since the different trajec- 
tories were not all recorded at the same frequency (96.5 
and 99.1 frames per second), we kept only the greater set 
(96.5 fps), containing 5 trajectories, and we split those 
into 526 short trajectories of 100 steps each. These tra- 
jectories arc non overlapping and one may view them as 
the result of 526 separate observations. Surprisingly, we 
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FIG. 6: Lipid granules diffusing in a yeast cell: log-log plot 
of the time-averaged second MME moment of the data from 
figure [5] as function of lag time (continuous lines), and Aot"'^ 
(dotted lines). Time is scaled in second, the ordinate is in 
/xm^. 



retrieve the exponent 0.41 ± 0.01 using the MSD, and 
the value 0.53 ± 0.02 from the second MME moment, as 
shown in figure [71 We repeated this analysis with a step 
size of 150 (350 trajectories) concluding that the choice 
of the step size 100 has no influence on the value of those 
coefficients. Since one of the trajectories (the magenta 
line in figures [5] and shows a much steeper slope, we 
excluded it for the rest of the analysis. 

An interesting observation is the following: assuming 
that the underlying stochastic process is indeed an FBM, 
relation (fTO)) for a = 0.41 predicts a value a' = 0.50 
for the MME statistics, in quite good agreement with 
the fitted value. This finding is quite suggestive in favor 
of FBM as the stochastic process governing the particle 
motion. 

Since the trajectories correspond to different granules, 
in different cells, we also studied them separately: each 
trajectory was split into stretches of 100 steps. For each 
granule, we plotted the regular and the MME ratios. 
They are somewhat noisy, but for each granule the MME 
ratio is clearly below the Brownian one (1.49): it ranges 
between 1.20 and 1.40. The regular moment ratio is 
slightly above the Brownian value (2), between 1.7 and 
2.5, as shown in figure 3 of the supplementary material. 
In the same figure we also plotted the ratio for the whole 
set of 100 steps pieces (thick lines), which give approxi- 
mately the same results as those obtained for individual 
trajectories. From these ratios, we obtain another clue 
pointing at an underlying FBM mechanism: the MME 
moment ratio is, on average, below the value for Brow- 
nian motion, and the regular moment ratio close to the 
Brownian value. These MME ratios are not very precise, 
but seem to range somewhat above the expected value 
for FBM with a = 0.41: equation (HI]) gives 1.21 ± 0.02. 
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FIG. 7: Lipid granules diffusing in a yeast cell: 526 sub- 
trajectories of 100 steps extracted from the experimental set 
of 5 trajectories, that are between 5, 515 and 19, 393 frames 
long. Ensemble averaged MSD (black o) fitted by a power 
law (q = 0.41, black line), and ensemble averaged MME (red 
o), fitted with a power law (a = 0.55, red line). We veri- 
fied that creating 350 trajectories of 150 steps instead of 100 
does not change the exponents obtained from the MSD or the 
second MME moment (x instead of o symbols). Since one 
of the trajectories had a steeper slope than the others, we 
repeated the same analysis without this trajectory, the new 
subset containing 445 trajectories of 100 steps, or 296 of 150 
steps (MSD in blue leading to a = 0.42, second MME mo- 
ment in magenta producing a — 0.51). Time is in seconds, 
the ordinate is measured in fj,m^ . 



The test with the growing sphere is, once again, some- 
what noisy, however, it clearly shows that the probability 
to be in a sphere, growing like t"^^, attains a constant 
value (see figure 4 of the supplementary material). This 
excludes the possibility that the process corresponds to 
diffusion on a fractal. 

The above analysis demonstrates that the tools pro- 
posed in this study allow us to classify the stochastic 
process underlying the motion of the measured single 
particle trajectories of the granules. We observe that this 
motion shares several distinct features with an FBM pro- 
cess. Namely FBM explains the finding of different scal- 
ing exponents of the MSD and the MME second moment, 
including their actual values connected by equation (jlOp . 
It is also consistent with a Brownian regular moment 
ratio, and an MME ratio lower than the Brownian one 
(compare figure 3 of the supplementary material). The 
recorded data were also shown to be incompatible with 
diffusion on a fractal. So what about CTRW as potential 
mechanism? The scatter between different single trajec- 
tories observed in the time averaged second moments is 
reminiscent of the weak ergodicity breaking for CTRW 
subdiffusion with diver ging characteristic waiting time, 
as studied in references [llf, [2^ . However an alternative 
explanation may simply be different environments and 
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granule sizes. It should be noted that even between suc- 
cessive recordings the cellular environment may change 
slightly, influencing the motion of the observed particle. 
The CTRW hypothesis however is not consistent with the 
moment ratio test: the expected ratio for a = 0.4 would 
be 3.38 for the regular one, and 2.50 for the MME, far 
above the observed values. 

Given the clues we obtained from the analysis, the ex- 
perimental data quite clearly point toward an FBM as 
underlying stochastic process. More extensive data ac- 
quisition is expected to allow more precise conclusions. 

V. CONCLUSIONS 

With modern tracking tools biophysical experiments 
provide us with the time series of single particle trajec- 
tories. Recently a growing number of cases have been 
reported in which the monitored particles exhibit subdif- 
fusion. An important example is the motion of biopoly- 
mers under cellular crowding conditions. While the mean 
squared displacement of these data, scaling like ~ , pro- 
vides the anomalous diffusion exponent a, the underlying 
physical mechanism causing this subdifFusion is presently 
unknown. As different mechanisms give rise to funda- 
mentally different physical behaviors influencing the par- 
ticle diffusion in a living cell, it is important to obtain 
information from experimental or simulation data other 
than the anomalous diffusion exponent, allowing us to pin 
down the specific stochastic process. We here introduced 
and studied several observables to analyze more quanti- 
tatively single particle trajectories of freely (sub)diffusing 
particles. For long trajectories with active motion events 
the latter may be singled out and our analysis performed 
on the passive parts of the trajectories [41|. As typical 
experimental data sets are relatively short, we here focus 
on the ensemble average obtained from a larger number 
of individual trajectories. The data were simulated on 
the basis of three subdiffusion models, these being con- 
tinuous time random walk with power law waiting time 
density, fractional Brownian motion, and diffusion on a 



fractal support. Moreover we analyzed two sets of exper- 
imental single particle tracking data, corresponding to a 
Brownian and a subdiffusive system. 

In particular we propose alternative measures to the 
usual fit to the mean squared displacement. Apart from 
obtaining the fourth order moment and construct the 
ratio (r^) / {r'^)'^ , these alternatives are: (i) mean maxi- 
mum excursion statistics that the particle has not trav- 
eled more than a preset distance up to time t. Its second 
and fourth moments, theoretically, scale with time the 
same way as the regular moments; however, they appear 
to reproduce more truthfully the actual subdiffusion ex- 
ponents. Constructing the ratio (?'max)/(''max)^ foi" these 
quantities provides additional information, that allows 
one to distinguish different subdiffusion mechanisms, (ii) 
The analysis using a growing sphere containing a certain 
portion of particles appears as a quite reliable method to 
obtain the (fractal) dimension of the underlying trajec- 
tory. 

An application to an experimental set proves the effi- 
ciency of those tests: the MME analysis is clearly more 
accurate than the classical MSD one, and with a modest 
data set we are able to collect several independent clues 
to identify FBM as mechanism to explain the motion of 
lipid granules under molecular crowding conditions. For 
long recorded time series the performance of the MME 
and regular moments analysis becomes comparable. 

From the discussion of simulations and experimental 
data is was shown that in order to understand the phys- 
ical mechanism of anomalous diffusion in a given set of 
data one needs to gather evidence from complementary 
measures, such as the ones proposed in this study. 
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